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Abstract 

The information on stellar parameters and on the stellar interior we can get by studying 
pulsating stars depends crucially on the available observational constraints: both seismic 
constraints (precision and number of detected modes, identification, nature of the modes) 
and "classical" observations (photospheric abundances, effective temperature, luminosity, 
surface gravity). We consider the case of (3 Cephei pulsators and, with the aim of estimating 
quantitatively how the available observational constraints determine the type and precision 
of our inferences, we set the stage for Hare&Hound exercises. In this contribution we present 
preliminary results for one simple case, where we assume as "observed" frequencies a subset 
of frequencies of a model and then evaluate a seismic merit function on a dense and extensive 
grid of models of B-type stars. We also compare the behaviour of x 2 surfaces obtained with 
and without mode identification. 

Session: Poster 



Tools 

In order to set the stage for Hare&Hound exercises, the following three main components 
need to be defined: 

• Theoretical predictions: The grid of models we use is betadat (Thirion &l Thoul 
2006). Stellar models and adiabatic frequencies are computed, respectively, with cles 
(Scuflaire et al. 2008a) and LOSC (Scuflaire et al. 2008b). The masses considered in 
the grid span the domain between 7.8 and 18.5 Mq, a metal mass fraction Z between 
0.01 and 0.025, an initial hydrogen mass fraction X = 0.70, and four values of the 
overshooting parameter (q; v) in the range 0-0.25. Frequencies of low-order oscillation 
modes of degree up to £ = 2 are computed for main-sequence models. 

• Observational constraints: we consider only seismic constraints, i.e. a subset of 
theoretical eigenfrequencies of a model M0 in the grid. The effects of rotation on the 
oscillation modes are not considered in this first step, thus all modes are assumed to 
be axisymmetric (m = 0). Concerning the degree of the observed modes, we consider 
the case where £ is unknown as well the one where £ is available as a constraint. 
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Table 1: Theoretical oscillation frequen- 
cies of model MO considered as observa- 
tional constraints, the uncertainty on the 
frequencies is assumed to be 0.1 ^Hz. 
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• Merit function: In order to compute a merit function at each point of the grid, we 
use a double optimisation procedure similar to the one extensively adopted in sdB aster- 
oseismology (see e.g. Brassard etal. 2001 and Charpinet et al. 2005). For each model 
in the grid we find the best global match between "observed" and theoretical frequen- 
cies by using a standard x 2 merit function. Then we study the properties of good-fit 
models looking at minima in the x 2 as a function of the stellar parameters/properties 
(e.g. location in an HR diagram, mass, central hydrogen mass fraction (X c ), mean 
density, . . . ). 



We consider as seismic constraints 4 oscillation modes of the model MO (see Table 1), with 
frequencies in the typical domain of the pulsation modes observed in (3 Cephei stars. The main 
parameters and properties defining MO are: M = 10 M , X c = 0.2, aov = 0, Z = 0.02, 
logTe ff = 4.34 and log L/L© = 4.02. 

We then compute the merit function on a sub-grid of models of betadat, where cnov ar| d 
Z are the same as in M0. The behaviour of x 2 for main-sequence models of different mass 
and evolutionary status is shown in Fig. 1. In the case where t is given as an observational 
constraint for all the modes (right panels), the properties and parameters of the input model 
M0 are easily recovered due to the appearance of an isolated global minimum in the x 2 
function. The increase of the number of x 2 minima in left panels (£ is unknown) compared 
to right panels allows to estimate the loss of information when mode identification is not 
available. Nevertheless, even if no mode identification is available, the frequencies of 4 
modes allow to constrain the parameter space in regions close to x 2 minima. 

The results of this simple test show that the approach presented here is viable tool to 
determine the number and precision of observational data required to constrain the properties 
of (3 Cephei pulsators. However, an extensive and thorough study is needed to investigate the 
effect of considering both additional uncertainties (e.g. on parameters such as overshooting 
and initial chemical composition, on rotational splittings and identification of m) as well 
as other constraints (e.g. luminosity, T e ff, \ogg, photospheric abundances, non-adiabatic 
constraints). 
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Figure 1: Colour-coded x 2 (logarithmic scale) in three different planes: log T e{{ -log L/Lq (upper panels), 
log T e ff-mean density (middle panels), and mass-X c (bottom panels). White dots represent the position of 
the input model MO in each parameter space. Right panels show % 2 computed including the identification 
of £ as a constraint, whereas in left panels i is considered unknown (i— 0,1 or 2). 



